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Abstract. More than 450 stars hosting planets are known today but only approximately 30 plane- 
tary systems were discovered around stars beyond the Main Sequence. The Penn State-Toruh Planet 
Search, putting an emphasis on extending studies of planetary system formation and evolution to 
intermediate-mass stars, is oriented towards the discoveries of substellar-mass companions to a large 
sample of evolved stars using high-precision radial velocity technique. We present the recent status 
of our survey and detailed characteristic for ^350 late type giant stars, i.e. the new results of radial 
velocity analysis and stellar fundamental parameters obtained with extensive spectroscopic method. 
Moreover, in the future we will make an attempt to perform the statistical study of our sample and 
searching the correlations between the existence of substellar objects and stellar atmospheric pa- 
rameters according to previous works which investigated the planetary companion impact on the 
evolution of the host stars. 
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INTRODUCTION 

The search for planets around giants represent an essential complement to "traditional" 
surveys, because they furnish information about properties of planetary systems around 
stars that are the descendants of the A-F main sequence (MS) stars with masses as high 
as ~ 5 M . As the stars evolve off the MS, their effective temperatures and rotation rates 
decrease to the point that their radial velocity (RV) variations can be measured with a 
few m s _1 precision. This offers an excellent opportunity to improve our understanding 
of the population of planets around stars that are significantly more massive than the 
Sun. Without such surveys it would be difficult to construct a broad, integrated picture 
of planet formation and evolution. 

Since 2001, about 30 such objects have been identified, including our five published 
detections with the Hobby-Eberly Telescope (HET, Niedzielski et al. 7, 4, 5). Our work 
has produced the tightest orbit of a planet orbiting a K-type giant identified so far (0.6 
AU) and the first convincing evidence for a multi-planet system around such a star [4]. 
Our recent discoveries [5] have also identified new multi-planet systems, including a 
very intriguing one of two brown dwarf-mass bodies orbiting a 2.8 M , K2 giant. This 
particular detection challenges the standard interpretation of the so-called brown dwarf 



desert known to exist in the case of solar-mass stars. Along with the discoveries supplied 
by other groups our work has substantially added to the emerging evidence that stellar 
mass positively correlates with masses of substellar companions, all the way from red 
dwarfs to intermediate-mass stars. 

In this paper we present the current status and forthcoming results from the Penn 
State-Toruh Planet Search (PTPS) around G-K giants with the HET since 2004. 

TARGETS SAMPLE 

High-precision radial velocity determination for MS stars more massive than ~ 1.3 M 
are not possible due to small number of spectral lines present in the spectra of these 
stars and their significant rotational broadening. However, it is possible to obtain precise 
RVs for intermediate-mass stars after they have left the MS and become subgiants or 
giants cool enough to develop spectra that are rich in narrow lines. Among the various 
planet detection methods in the case of stars more massive than ~ 1 M Q RV searches 
provide the most efficient way to do it, while for stars more massive than ~ 2 M they 
are practically the only way. 

The whole PTPS targets sample of nearly 1000 late-type stars observed with this 
survey (Fig. 1) comprises about 350 giants from the giant clump, about 350 giants and 
subgiants, and roughly 250 evolved stars in the upper envelope of the MS (aging dwarfs). 
In terms of using the RV technique, this star sample is limited to 0.4 m < (B — V ) < 1 .2 m 
and concentrated around 1.5 m above the MS. Moreover, all PTPS targets are uniformly 
distributed in right ascension and declination on the sky. The magnitude upper limit of 
the stars in the V band is 1 l m and usually (> 66 %) the program stars are fainter than 8 m . 

OBSERVATIONS AND DATA REDUCTION 

Observations for the PTPS have been obtained with the 9.2 m HET [8] equipped with the 
High Resolution Spectrograph (HRS, R=60000, Tull 14) in the queue scheduled mode 
[9]. The spectrograph was fed with a 2 arcsec fiber. The spectra consisted of 46 echelle 
orders recorded on the "blue" CCD chip (407-592 nm) and 24 orders on the "red" chip 
(602-784 nm). A typical signal-to-noise ratio was 200-250 per resolution element. The 
basic data reduction was performed in a standard manner using the IRAF 1 tasks and 
scripts. 

The observing strategy was focused on long-term variations in RVs which were 
measured using the standard I2 cell calibration technique [2]. Details of our survey, the 
observing procedure and data analysis, have been described elsewhere [7, 6]. 



IRAF is distributed by the National Optical Astronomy Observatories and operated by the Association 
of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 
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FIGURE 1. The PTPS survey stars: clump giants (circles) contains stars of various masses over a range 
of evolutionary tracks; intermediate-mass dwarfs (squares) which are about to leave the MS; giants and 
subgiants (triangles). Background stars (grey points) are from The Tycho-2 Catalogue and are not included 
in the PTPS. 

STELLAR CHARACTERISTICS AND RESULTS 

Late type giants are known to exhibit various types of intrinsic variability induced by 
pulsations and/or surface inhomogeneity connected with rotation. These effects have to 
be ruled out before substellar companion interpretation can be considered. Therefore, 
proper interpretation of results obtained from RV survey of G-K giants requires detailed 
knowledge of their physical parameters. 

Atmospheric parameter determination 

An atmospheric model of the star is characterized by effective temperature T e ff, 
surface gravity logg, microturbulence velocity v t , and metallicity [Fe/H]. To obtain 
these parameters from high resolution spectra, we used a purely spectroscopic method 
developed by Takeda et al. [13, 1 1, 12] which is based on the measurements of equivalent 
widths (EWs) of iron lines and resulting from the assumption of LTE. We were able 
to find and measure the EWs of ~300 Fe I and Fe II lines in our spectra. To avoid a 
perturbation of the mean iron abundance (and thus [Fe/H]) dependency on the excitation 
potential and EWs, we selected the lines not stronger than 150-200 mA. 

Figure 2 presents the results determined from the method mentioned above and 
the dependencies between atmospheric parameters (r e ff, log g, v t , and [Fe/H]) in the 
case of clump giants. The consistency of the solutions is clearly visible and remains 
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FIGURE 2. Dependencies between four atmospheric parameters for PTPS clump giants: 
(a) r e ff versus logg, (b) T e ff versus [Fe/H], (c) logg versus [Fe/H], and (d) logg versus Vt. 
Grey points indicate eight stars with substellar-mass companions discovered by the PTPS. 



in agreement with common features of Red Giant Clump (RGC) stars. The average 
uncertainties of the following parameters are oT e ff = 14K, dog (g[cms~ 2 ]) = 0.05, 
ovi = 0.08kms _1 , and a [Fe/H] = 0.07. Only for individual stars these values are slightly 
higher. 



Integral parameters determination 

The knowledge of T ef f, log g, and [Fe/H] together with Hipparcos/Tycho data allows 
the estimation of stellar bolometric magnitudes Mboi and luminosities log(L/L ). The 
rest of necessary quantities, such as intrinsic colour (B — V)o and bolometric correction 
BC, can be determined for each star using empirical calibration taken from Alonso et al. 
[1]. To obtain the stellar mass, radius, and age, i.e. the parameters constraining the stellar 
evolution, the position of the each single star should be compared with appropriate 
evolutionary tracks and isochrones on the H-R diagram [3, 10], although the accurate 
estimation of these parameters is made very difficult due to the complex nature of the 
RGC region in the H-R diagram. 

Figures 3 and 4 show the distributions of masses and radii estimated from the gen- 
eral definitions based on T e f f , logg, and log(L/L ). It is worthwhile noting that these 
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FIGURE 3. Masses distribution of the PTPS clump giants. Most of the stars fall in the 1-3 M range. 
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FIGURE 4. Radii distribution of the PTPS clump giants. Most of the stars fall in the 8-13 R range. 
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estimates do not take into account the stellar composition, which is a very important fac- 
tor influencing the stellar evolution. Nevertheless, both distributions under consideration 
confirm that we study the objects in the range 1-3 M and 8-13 R©. 

The typical systematical errors are clog(L/L ) = 0.4, oM = 0.3 M , and oR = 
0.8 R . However, we should mention that for many cases the uncertainties are much 
higher because of the moderate quality of stellar photometry available in the astronomi- 
cal databases and/or high uncertainties in Hipparcos parallaxes. 



SUMMARY AND CONCLUSIONS 



Since 2004, conclusive multi-epoch observations have been gathered for about 800 stars. 
The large number of data allow for a preliminary assessment of possible companion 
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FIGURE 5. The first ten substellar-mass companions discovered by the PTPS in function of the distance 
from the hosting stars. The symbol sizes are proportional to planet mass. 



candidates, and, depending on the criteria adopted, the substellar-mass companion fre- 
quency may reach up to 30% in the case of the whole PTPS sample. 

Until today seven stars from PTPS survey have been confirmed to be orbited by 
planetary systems. The first ten substellar-mass companions discovered by the PTPS 
are shown in Fig. 5. In addition to the published planet and brown dwarf companions, 
we also present two new planets around BD+20 274 and HD 240237 (Gettel et al. in 
preparation), and one around HD 17028 (Nowak et al. in preparation). 

The precise knowledge of stellar parameters is also important in terms of the statistical 
correlations between planet occurrence frequency, their orbital parameters, and stellar 
mass or metallicity. The forthcoming study (Zielihski et al. in preparation) shall give a 
new set of data, and together with previous work [15, 16], it will bring more information 
for statistical analysis and improve the general knowledge concerning the evolution of 
stars and planetary system formation scenarios. 
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